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Ekpyrotic bouncing cosmologies have been proposed as alternatives to inflation. In these scenarios, 
the universe is smoothed and flattened during a period of slow contraction preceding the bounce 
while quantum fluctuations generate nearly scale-invariant super-horizon perturbations that seed 
structure in the post-bounce universe. An analysis by Tolley and Wesley (2007) showed that, for 
a wide range of ekpyrotic models, generating a scale-invariant spectrum of adiabatic or entropic 
fluctuations is only possible if the cosmological background is unstable, in which case the scenario 
is highly sensitive to initial conditions. In this paper, we analyze an important counterexample: a 
simple action that generates a Gaussian, scale-invariant spectrum of entropic perturbations during 
ekpyrotic contraction without requiring fine-tuned initial conditions. Based on this example, we 
discuss some generalizations. 
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I. INTRODUCTION 

Observations of the cosmic microwave background 
from the Wilkinson Microwave Anisotropy Probe 
(WMAP) [T], the Planck satellite H 13], the Atacama 
Cosmology Telescope (ACT) ^ and other experiments 
have shown that the primordial scalar (density) fluctua¬ 
tion spectrum is adiabatic and nearly scale-invariant with 
nearly Gaussian statistics. Inflation m has been sug¬ 
gested as a mechanism for generating perturbations with 
these properties, though, it is known that to do so it 
requires rare initial conditions [El IS] and results in a mul- 
tiverse of outcomes [Mill- 

Bouncing cosmologies with a period of ultra-slow 
(ekpyrotic) contraction have been proposed as alterna¬ 
tives. In these theories, smoothing contraction occurs 
because the energy density of a scalar field with equation- 
of-state e > 3 (where e = 3(p -I- p)/p with p being 
the pressure and p the energy density) grows to domi¬ 
nate all other forms of energy, including inhomogeneities, 
anisotropy and spatial curvature m- A key advantage 
compared to inflation is that the ekpyrotic mechanism 
does not lead to a multiverse. 

The currently best understood way to produce den¬ 
sity fluctuations in the ekpyrotic theory involves two 
scalar fields that generate a scale-invariant spectrum of 
entropy perturbations. After the ekpyrotic smoothing 
phase, these perturbations convert into a scale-invariant 
spectrum of adiabatic perturbations [THlllflj . In the first 


examples discussed in the literature, the background cos¬ 
mological solution describing the evolution of the two 
fields along the potential energy surface is unstable, 
which means finely-tuned initial conditions are required 
to begin the ekpyrotic phase [20H23]. Tolley and Wesley 
[24] analyzed the dynamics of a wide class of contracting 
cosmological models that generate scale-invariant adia¬ 
batic or entropic perturbations and suggested that the 
problem may be generic. More specifically, they showed 
that the cosmological background solutions are not at¬ 
tracted to a fixed point and, from this, concluded that 
the models are highly sensitive to initial conditions. 

In this paper, we present simple ekpyrotic models that 
generate a scale-invariant, nearly Gaussian spectrum of 
density perturbations but do not require fine-tuning of 
initial conditions. Although these models belong in the 
class considered by Tolley and Wesley, we show that the 
background solutions are attracted to a fixed-curve along 
which scale-invariant fluctuations are generated. The ex¬ 
istence of such a fixed-curve is sufficient to ensure that 
the observational predictions are insensitive to the choice 
of initial conditions. In other words, being attracted to 
a fixed-point is not necessary to avoid fine-tuning. 

In Section [n] we summarize the general argument that 
suggests the need for finely-tuned initial conditions. In 
SectionjnTj we review a simple ekpyrotic model for which 
we find a fixed-curved but no fixed-point attractor. In 
Section |IV[ we describe how to construct more general 
examples that also avoid the need for fine-tuning of ini¬ 
tial conditions. Finally, we discuss the implications for 
cosmology. 


aijjas@princeton.edu 



2 


II. SCALING SOLUTIONS, 
SCALE-INVARIANCE AND INSTABILITY 


Scaling solutions are solutions to the equations of mo¬ 
tion for which there exists a set of field variables such 
that all contributions (treating the kinetic and poten¬ 
tial energy densities as distinct) to the total energy den¬ 
sity scale identically with time, keeping their fractional 
contributions constant. Scaling background solutions are 
particularly important because they are exactly solvable 
and can yield a scale-invariant spectrum of perturbations. 

In Ref. [23] , Tolley and Wesley presented an instability 
argument that applies to contracting, scaling solutions 
derived from two-derivative, two-field actions 

S = j d^Xy/^g^R (1) 

that possess a continuous symmetry generated by a pa¬ 
rameter K such that 

d^ 

^=ea($), S^e-S. (2) 

Here is the spacetime metric with (—h -l--|-) signa¬ 
ture convention, R is the Ricci scalar, ^a($) is a function 
of the two scalar fields = {$a} where a = 1,2, 
is the metric on field space, and H($) is the potential 
energy density; reduced Planck units (SttGtv = 1 where 
Gat is Newton’s gravitational constant) are used through¬ 
out. This symmetry guarantees the existence of a set of 
field variables ($ 1 , $ 2 ) t o’) such that the Lagrangian 
density can be rewritten as 

C=Ir- lidaf - Ifiamr - Hoe-“^h(a), (3) 


where c is a real constant and Vb < 0 (that is, the 
ekpyrotic potential V{(j),a) is negative). Along the back¬ 
ground solution, /(O) = h(0) = 1 and cr = 0. 

It proves useful to introduce the dynamical variables 


{w,x,y,z) 


VQR ’ V6H 


a\/-Voh{a)e 

Vln 



where r < 0 is conformal time running from large nega¬ 
tive to small negative values, a prime denotes a derivative 
with respect to conformal time, a is the scale factor and 
R = a'/a is the conformal Hubble parameter. The four 
variables are dimensionless using reduced Planck units. 
With these variables, the Friedmann-Robertson-Walker 
(FRW) equations of motion become 

1> (" - 7m) “ vf 

1> (*++vf lb 

Z,N = VQx. (7) 



FIG. 1. A schematic trajectory of a solution through field- 
space: The field variables {4>,cr) were constructed so that, 
along the field-space trajectory in Eq. ( |10[ ), only 0 varies and 
CT = 0. Perturbations Scj) are tangent to the trajectory and are 
adiabatic; perturbations Sa are orthogonal to the trajectory 
and are entropic. 


Here we introduced the dimensionless time variable 
N = In a that denotes the number of e-folds of ekpy¬ 
rotic contraction and runs from large positive to small 
positive values. We eliminated y using the Friedmann 
constraint 


2,2 2 1 
w + X — y =1. 


( 8 ) 


The equation of state takes the simple form 

W 


e = l- ^ + x^). 


(9) 


From the Friedmann constraint, we also see that the 
square of each variable w, x, and y is a fractional con¬ 
tribution to the total energy density: is the (()-kinetic 

energy; is the cr-kinetic energy; and y^ is the potential 
energy. The equation-of-state parameter e is the sum of 
kinetic energies. 

As can be verified by direct substitution, the back¬ 
ground Eqs. (5][7) admit a fixed-point solution. 


{w,x,y,z)= (^, 0 , 


(10) 


in which a — 0 with and only (f> is changing, provided the 
constraint (0) = (0)/(c^ ~ 6). Obviously, this 

is a scaling solution since any fractional contribution to 
the total energy density, and is constant. 

The cosmological background solutions correspond to a 
field-space trajectory like the one shown in Fig.j^ Pertur¬ 
bations of this trajectory can be decomposed into those 
along the red curve (adiabatic perturbations) and per¬ 
pendicular to it (entropic perturbations). 

The instability argument connects the stability of the 
solution in Eq. ( [l0| with the spectral indices of its per¬ 
turbations. The basic idea makes use of the fact that, 
since z = a = 0 in the background solution, the second 
order action and, hence, the perturbation spectra derived 
from it are determined by a few parameters 


( 11 ) 
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FIG. 2. The instability of the scaling solution correspond¬ 
ing to the Lagrangian density in Eq. (12l: The background 
cosmological solution corresponds to a trajectory along the 
ridge of the potential as indicated by the red arrow and red 
curve. Quantum fluctuations along the trajectory produce 
a blue spectrum of adiabatic perturbations and fluctuations 
normal to the trajectory produce a scale-invariant spectrum 
of entropic perturbations. After the ekpyrotic phase, the en- 
tropic perturbations convert into adiabatic perturbations due 
to the bending of the trajectory curve (not shown here). The 
fact that the background trajectory is a ridge means that it is 
unstable if the initial conditions are sufficiently far from the 
ridge. 


Linearized around the fixed-point in Eq. (101, the back¬ 
ground Eqs. (5 [ 7 1 reduce to a matrix equation, where the 
parameters in Eq. (1^ determine the eigenvalues of the 
matrix. Tolley and Wesley’s analysis showed that any 
combination of the parameters in Eq. (Ill that results in 
a scale-invariant spectrum of perturbations (entropic or 
adiabatic) renders the background solution dynamically 
unstable to perturbations in the sense that the matrix 
associated with the linearized system has at least one 
negative eigenvalue. In a contracting universe, a nega¬ 
tive eigenvalue means a dynamically unstable direction 
in (w-a:-z)-space. 

An example is given by the Lagrangian density 


( 12 ) 

where ci,C 2 are positive-definite constants and Vq < 0 . 
The corresponding FRW equations of motion admit a 
scaling solution with = Ai\n\T\ + Bi and ciAi = 
C 2 A 2 that has been shown to generate a scale-invariant 
spectrum of entropic perturbations [HI [11]. According 
to the instability argument, it should have an unstable 
background which in this simple case {Gab = Sab) can be 
depicted as in Fig. 

For the purpose of illustration, we briefly outline how 
the instability emerges from a negative eigenvalue of the 
linearized system. The change of variables. 


_ C2$l -b Cl$2 
\/c? + 

Ci$i - C2$2 , 

+ O-Q 


(13) 

(14) 


with (To = 2 ln(c 2 /ci)/(cf -b c^), brings the Lagrangian 
density in Eq. (HI to the form of Eq. (|^ . The coupling 
function to the kinetic energy of 4>, f{a), and the coupling 
function to the potential energy of (f>, h{a), are given by 

/(^) = 1, (15) 

h{a) = l+^a'^+ 0{a^), (16) 

and the parameters, c and Vq are defined such that 



Linearizing the background Eqs. (^) about the fixed- 
point in Eq. (101, the perturbations (Sw, i5a;, i5z) = {w — 
-^,x,z) satisfy 


M = 


V 


Sw,N 

\ 5w 

6x,n 

1 = M • 1 5x 

Sz,N 

) 

3 

-6) 

0 

1 

2 



76 


(19) 


• ( 20 ) 


0 


The eigenvalues of M are (c^ — 6) /2 and 
((c^ — 6) ± \/9c^ — GOc^ -b 36) /4. Note that the 

smallest eigenvalue is negative for ekpyrosis, as is clear 
from substituting Eq. ( [Toj ) into the equation of state, 
Eq. (§: e > 3 requires c > y/E. Therefore, as the 
universe contracts, N decreases, and perturbations 
along the eigenvector corresponding to the negative 
eigenvalue grow so that the system is carried away from 
the fixed-point solution in Eq. (10). In this case, the 
negative eigenvalue means that the initial conditions for 
the fields must be fine-tuned to lie close to the trajectory 
or else the fields will evolve far-off course as illustrated 
in Fig. [§ 


III. EKPYROSIS AND SCALE-INVARIANCE 
WITHOUT FINE TUNING 


In this section, we describe the case where the nega¬ 
tive eigenvalue exists but is physically irrelevant. As we 
will show, the occurrence of the negative eigenvalue only 
means that the attractor is a fixed-curve rather than a 
fixed point. 

We consider the Lagrangian density first discussed by 
Li in Ref. [25], 


C = \R-l{d^)^- 


a = 


\/4 +4 


le-^^{dx)^-Voe-^^ 


( 21 ) 
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FIG. 3. The Lagrangian density in Eq. (211 with V = 
0 and A = \/6: Green curves and red arrows map onto 
each other according to the transformations in Eqs. ( |38[ 
|41[ |. Solutions follow the blue arrows. FIG. 3.a. The 
trajectories in (ip, x) variables: For y = Q, (w, x, z) = 

(^ip'/ix^U), e~^*(x-Xo)) ■ The green 


curve, {w,x,z) = (1,0,2), is the fixed-curve attractor. The 
red arrow points to the fixed-point (w,x,z) = (1,0,0). 
FIG. 3.b. The trajectories in (0, a) variables: For y = 


0, {w,x,z) = f((y)(t>'/(V&'H), a'a-'j . The 


green curve, (w,x,z) = (sech(A «/2), — tanh(A z/2), z), is the 
fixed-curve attractor. The red arrow points to the fixed-point 
(w,x,z) = (1,0,0). 


where A is a positive and Vq < 0. The model involves 
an ekpyrotic field, ip, with a negative potential, similar 
to ordinary, single-field ekpyrosis. The novel feature is 
the non-canonical, exponential coupling to the massless 
spectator field, y. We begin with the simple case where 
V{ip) = 0. This corresponds to the borderline ekpyrotic 
equation of state e = 3. Then, we generalize to V{ip) ^ 0 
(e > 3) and provide a full, analytic treatment. 


A. V(ip) = 0 


The basic idea is captured in Fig. [^, which illustrates 

background trajectories corresponding to different initial 

0 

conditions in the space (w,x,^,z) = 


JL 


e-2^x' 


Vsn 


,e ^'^(x-Xo) 


( 22 ) 

We use the superscript ^ to distinguish quantities ex¬ 
pressed in (ip, x) variables from those expressed in the 


special ((p, a) variables that were used to derive the orig¬ 
inal instability argument. 

Any set of initial conditions for a, ip, ip', x, and x' cor¬ 
responds to a particular point on the cylinder. The back¬ 
ground solution follows the blue arrow originating at this 
point. There are two special initial conditions at the 
points (w, X, z) = (±1,0,0); the one with the -I- sign cor¬ 
responds to the red arrow in Fig. [^. These two points 
are special because the blue arrows vanish here, i.e., if 
the background solution starts here, it stays here. Hence, 
these are fixed-point solutions. 

For any other initial conditions, (i.e., any other point 
on the cylinder), the blue arrows carry the background 
solution toward the green curve. The green curve is 
therefore a strong attractor for generic initial conditions. 
Henceforth, we call such an attractor curve a “fixed-curve 
attractor.” Solutions evolving along this curve are scaling 
solutions since w,x, and y are constant. They describe a 
universe dominated by the kinetic energy of the adiabatic 
field. Ip] the entropic fielH y, is constant. For the case 
(A = V&) depicted in Fig.|^ these solutions can be shown 
to generate a scale-invariant spectrum of perturbations in 
the entropic field, y. 

As for the negative eigenvalue associated with the lin¬ 
earized equations of motion, in the cases discussed in this 
paper, it only indicates the existence of a fixed-curve at¬ 
tractor instead of a fixed-point attractor. The existence 
of a fixed-curve attractor means there is no need for fine- 
tuning of initial conditions. 

One can describe the same dynamics in the special 
((p, cr) variables (see Fig. Since the Lagrangian den¬ 
sity in Eq. © has the shift symmetry of Eq. ([^ as¬ 
sumed by the instability argument, it can be put into 
the form of Eq. (§ through a variable transformation 


iP 





y = — e 2 tanh 

A 



+ V’o, (23) 

+ X0, (24) 


where ipo is a real constant. Substituting this transfor¬ 
mation into Eq. (21) yields the Lagrangian density in 
Eq. ([^ with /(cr) = h(a) = cosh^(Acr/2) and c = A. 
The two Lagrangian densities describe the same theory 
in different field variables: the cylinder on the right is a 
twisted version of the one on the left. It is clear from 
the blue arrows that the green curve in Eig. is an at¬ 
tractor, just like the green curve in Eig. Solutions 
along the green-curve attractor generate a scale-invariant 
spectrum of entropic perturbations. 


B. V(ip) ^ 0 


We extend our analysis to the more general form of 
the Lagrangian density in Eq. (21) that applies both to 
V(ip) = 0 and to V(ip) ^ 0. 
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Extremizing the action with respect to variations of 
the fields {ip,x) yields the field equations 


iIj” + - XVoe-^'f’a^ + ^ g, 

X" + 2Hx' - AV-'x' = 0. 

The Friedmann constraint is 


(25) 

(26) 

(27) 


Using the variables defined in Eq. (22), these equations 


can be recast as the autonomous, dynamical system 


= 3(u}^ - 1) ( w - ^ ) - 


W,N 


X,IS! = 3x(w"‘ + x"‘ — 1) + \/ -Xwx, 




Z,N = —\I-Xwz + V6x. 


(28) 

(29) 

(30) 


If U 7 ^ 0, this system admits three fixed-point solutions 
at {w,x,y,z) = 


(- 1 , 0 , 0 , 0 ) 
(+ 1 , 0 , 0 , 0 ) 




(31) 

(32) 

(33) 


all of which are unstable {i.e., associated with a nega¬ 
tive eigenvalue). The third fixed-point solution given by 
Eq. (p^ bisects two fixed-CMr?;e solutions 


{w,x,y,z) = ( 


- 1,±Z 


(34) 


with Z (X that generate a scale-invariant spectrum 

of entropic perturbations, as shown in Ref. [25] . 

If U = y = 0, X^ must be 6 in order for the fixed-point 
in Eq. ([M|) to be a solution. For A = -x/O, this coincides 

and corresponds to the 


with the fixed-point in Eq. 
red arrow shown in Fig. Eq. ( [M| parameterizes the 
vertical, green, fixed-curve attractor. 

Changing variables {-ip, x) t o') as defined in 
Eqs. (23) and (24), the Lagrangian density in Eq. (21) 
takes the form of Eq. ([^ with f{a) = h{a) = 
cosh^ {Xa/2) ,c = A, and Vq = Repeating the 

same analysis in the new variables defined in Eq. ([4]), the 
equations of motion become 

WjM = + x"^ — 1) (w -^sech (Az) 

V V6 

— ■\/6ctanh(Az) a;in, (35) 

x,N = 3(w^ + — 1) + ^/^Atanh(Az)^ 

+ •\/6ctanh(Az )w^, (36) 

z,N = VGx. (37) 


The variable transformations as defined in Eqs. (23) 


and (24) imply the following relationship between the 


variables {w,x,y,z) i—>■ {w,x,y, z): 


w = sech ( —w — tanh I 
X = tanh ( —zj w + sech ^ — z ) x, 


y = y, 


(38) 

(39) 

(40) 

(41) 


If U = 0, these transformations quantify how to “twist” 

Fig.jl. to generate Fig. §1. _ 

The fixed-point solutions in Eqs. ( M]|M ) are given in 
the new variables {w, x, y, z) as 


(- 1 , 0 , 0 , 0 ) 
(+ 1 , 0 , 0 , 0 ) 




(42) 

(43) 

(44) 


and the fixed-curve solutions in Eq. (34) are 





with Z = (2/A)sinh ^ (^XZ/ 2 '^. 


(45) 

These fixed-curves lie 


on the surface of the cylinder +x^ = A^/6. For U = 0 
and A = \/6, Eq. (45) corresponds to the twisted green 


curve that is confined to the surface of the unit cylinder 
Fig. |3 }d. 


Direct substitution verifies that the curves in Eq. (p34]) 
and Eq. 


are solutions to the background equations 
given in Eqs. ( 28]|30 ) and Eqs. ( 35]|37 ) for both U = 0 
with A = \/6 and for V ^ 0. 

To show the existence of a negative eigenvalue, we lin¬ 
earize the equations of motion about the fixed-points. 


in Eq. (33) and (44), respectively, for the two sets 
of variables. Linearizing Eqs. ( |28||30 ) about Eq. (33) 
yields a matrix equation like that given in Eq. (19) with 
{Sw, 6x, 6z) = (w — A/-\/6, x,z) and 


M = 


V -3 0 

0 A2-3 


V6 -4^ 


(46) 


Similarly, linearizing Eqs. ( [36||37| about Eq. ( [44| ) yields 
Eq. (19) with {6w, Sx, Sz) = {w — A/-\/6, x, z) and 


M = 



(47) 


Both M and M have eigenvalues 
{—A^/2, A^/2 — 3,A^ — 3 }, the first of which is negative. 
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In the first set of variables, the eigenvector corresponding 
to the eigenvalue — A^/2 is parallel to the unit vector 
in the z-direction, which is tangent to the fixed-curve 
solution in Eq. (34). In the second set of variables, 
the eigenvector associated with the eigenvalue — A^/2 
is parallel to a linear combination of unit vectors x, i, 
namely z — X^j(2^/6) x, that is tangent to the fixed-curve 
solution in Eq. ( [4^ . 

For the case, E = 0 and A = y/6, these eigenvectors 
are tangent to the green fixed-curves in Fig. at the red 
arrows. The existence of a negative eigenvalue in this 
model is harmless, since it only means that the system 
is attracted to a fixed-curve solution (instead of a fixed- 
point solution) that generates a scale-invariant spectrum 
of entropic perturbations. 


C. Further Generalizations 


Although the remainder of this work will consider ac¬ 
tions with the shift symmetry in Eq. ([^, our results can 
be generalized to cases without shift symmetry. For ex¬ 
ample, the field contribution to the shift-symmetric La- 
grangian density in Eq. (21) is a special case of the more 
general Lagrangian density 


£=^R-^{di}f-^e ^^{dx)^-{l + r{xj)Voe ^'^+q{x) 

(48) 

with fi = X and r(x) = q{x) = 0. The addition of q{x) 
and r{x) breaks the shift symmetry since 

E = (l + r(x))Eoe-^’^ + g(x) 

—>■ ^1 -I- r{e^'^x)'j Voe~^^~'^ + r{e^'^x) 

+ e"”E. (49) 


If ^ = A, the ekpyrotic Lagrangian density in Eq. (48) 
admits a scaling solution that is a fixed-curve attractor 
with = 0 and that generates a scale-invariant spec¬ 
trum of entropic perturbations. This is due to the fact 
that, as x! 0i ^(x) ^md q{x) approach constants 'r(xo) 
and q{xo)- The first, r(xo), can be reabsorbed into Vb, 
and the second, q{xo), is negligible along the fixed-curved 
attractor. 


IV. CONSTRUCTING NEW MODELS 

In this section, we derive the most general ekpyrotic, 
two-field Lagrangian density with shift symmetry that 
admits scaling solutions which are either fixed-point or 
fixed-curved attractors and generate a scale-invariant 
spectrum of entropy perturbations. 

First, we consider the Lagrangian density in Eq. ([^ 
with arbitrary parameters and couplings, 

{Eo < 0,cG K,/i(a) > 0,/(ct) > 0}. (50) 


In the Appendix, we show that the combined conditions 
of shift symmetry, scaling solution, fixed-curved attractor 
and scale-invariant spectrum of perturbations imply the 
following properties: 

PI: lim /(cr) = oo monotonically; 

|(t|—>- oo 

P2: lim h{a) oc 

|(7|—>-00 


P3: lim {w,x,y,z) 

|f7|—>-00 



l,-sgn(^)oo ; 


P4: \fj.\ > Ve. 

Property PI says that the coupling f{a) must grow with¬ 
out bound. Property P2 constrains the form of the po¬ 
tential energy density to be exponential at late times. 
Property P3 defines the scaling solution. It implies w = 0 
so that the (f) field is fixed; furthermore, since the back¬ 
ground equations in Eqs. (§71 depend explicitly on a 
which will vary, it also implies that the scaling solution 
is a fixed-curve (rather than a fixed-point) in (w-x-z)- 
space. Property P4 is necessary for ekpyrosis (e > 3) as 
follows from substituting this solution into the equation 
of state, Eq. ([^. 

The example from the last section has these four 
properties. At late times, the fixed-curve attractor in 
Eq. ( [4^ goes to z = cr —>■ ±oo. In this limit, /(cr) = 
h{a) = cosh^(Acr/2) is dominated by the single exponen¬ 
tial /4. For example, if A > 0 and a —>■ —oo, choosing 
/i = A in the solution in property P3 reproduces Eq. (45) 


at late times. Similar arguments apply for the different 
combinations of sgn(A) and sgn(cr). 

Assuming properties PI thru P4, the only remaining 
degrees of freedom are the parameters, c, Eq, and the late¬ 
time behavior of /(cr), modulo property PI. We show 
now that, given these four properties, it is possible to 
obtain a scale-invariant spectrum for the entropic modes 
but not for the adiabatic modes. 

We perturb Einstein’s equations about the fixed-curve 
solution specified by Property P3, working in the longi¬ 
tudinal gauge [HHZ] where the metric takes the form 


ds^ = a2 24>)dr2 + (1 - 2^)dx^) . 


(51) 


Since </>' = 0 along the background solution in prop¬ 
erty P3, the quantity Qs = \/is automatically 
gauge-invariant and represents the entropy perturbation; 
the Mukhanov-Sasaki variable Qa = 6a + {a'/%) $ is 
also gauge-invariant and represents the adiabatic pertur¬ 
bation [551 HH] ■ Property P3 implies that the equation 
of state e = y?/2 and, therefore, the conformal Hubble 
parameter is = (1 — e)(—r) < 0. Then, the mode 
functions = a and Ug = aQs can be shown to 
satisfy 



(-r)2 


(-r)2 


1l(j 

Us 


(x [3 

'y 6 

- - -Ua + 7 - zu'^, (53) 

(-r )2 (-t) 
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where the the background-dependent quantities can be 
derived, for example, from the expressions in Ref. [30]: 


r = -2 


= -2 


- 4 


- 6)' 


(^^2 - 2)2 _ 2)2 f 


- 4 

(^2 _ 2)2 + 


^2-6 - 6) 
(^ 2 - 2)2 / 

f jaa 


- 2 ) / 

2 o .,2 


a = 4' 
/3 = 2 
7 = - 

<5 = -/ 


1 

to 

1 

2)2 

V 7 J 

1 

6 

c 

m(m^- 

-2)2 

Vf ' 

o 

6 

c 


-2)- 



- 6 

c 


2/i 

^ (^2 _ 2)2 ^ 

, - 6 c /,„ 

‘(m^-2)V7 / 


/r (/r2 - 2)^ yj’ 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 


These variables are generally time-dependent because 
/ = /(cr) is a function of r. Substituting the emres- 
sion for into the definition of x given in Eq. (Wl), we 
find 




(60) 


Using this expression, Eqs. ( |54]|5^ can be rewritten as: 
- 4 (^^-6)^ 


r = -2 

r = -2 
1 

~4 
a = 4 


^ 2(^2 - 2)2 / 
— Q 


(^2 - 2)2 

/r2 — 4 

(m2-2)2^%2(^2_2)7 

j(-r)j +5j(--). 


/' 


/ 


/? = 2 


— 6 c — 6 c f 

^(^2 _ 2)2 yj + ^(^2 _ 2 ) 77 / 

y? — & c 


{-t), 


^i{p^-2) 77’ 


7 = —4 


— 6 


^ (^2 _ 2)^ 77’ 


(5 =-7 


(61) 

(-t) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


Eqs. (521 and (53) are a coupled, linear system of dif¬ 
ferential equations which must be solved as r —)■ 0“ to 
find the spectra. Depending on the growth rate of the 
coupling function f, different terms in Eqs. ( 6T||66 ) come 
to dominate in this regime. For example, the first term 
in 9'^ always dominates over the second since 7 —t oo as 
T —>■ 0. By contrast, the relative sizes of the different 
terms in 0® depend on the magnitude of f' /f. 

For clarity, we define the symbol to mean 


I \i ^ In/ M T |A(t)| d\n\B{T)\ 


> 0 . 

(67) 


Similarly, we define to mean 

( 68 ) 

For example, l/(—r)^ > l/(—r) and 3/(—r) ~ 2 /(—t). 

There are three qualitatively different cases to con¬ 
sider: fast growth, |/,„ //| > 1, slow growth, |7,„ /f \ < 1, 
and “just-so” growth, \f,a/f \ ~ 1- 


A. Fast growth: |/,ct//| ^ 1 

If 17:0-/7I ^ 1> then I/V7I ^ l/(-'r)- Let us assume 
that 


7?-<3 7-7 < S 


(-r)2’ i-r) 


(-r)2 


(69) 


so that to leading order, the entropic mode evolves inde¬ 
pendently: 


u:_u£, ,...0. 


2 V / 2 V/, , 

For the mode function we find the solution 

df 


= \/7M ci(fc) 


-l/k 


7('r) 


C2(fc) 


(70) 


(71) 


where Ci(fc) and C 2 (fc) are constants of integration. 
Choosing ci{k) and C 2 (fc) so that Us and u(, match the 
Bunch-Davies solution, (l/-\/2fc)e“®^'^, at horizon cross¬ 
ing, —T = 1/fc, yields 


ci{k) = 
C2{k) = 


e- (f (-l/fe) + 2»fc7(-l/fc)) 

2y/2kf{-l/k) 


72 A:/(-l/fc)' 


(72) 

(73) 


For fast-growing f, at late times the integral in Eq. 0 
is very closely approximated by (l/fc)(l/7(—1/fc)). With 
Eqs. (72) and (73), the entropic mode function is given 
by 


Us = E(fc)v7(7) 


(74) 


where 


. _ U (-7 (-l/fc) + (2-2z)fc/(-l/fc)) 

2y2fc3/2(/(-l/fc))3/2 ■ 

Substituting this result into the right side of the adiabatic 
equation, Eq. (52), we find that satisfies 


u'' + 2 _ —u - 

(/r2-2)^ (-r)2“'^-%(^2_2)2 (-r)2 


(76) 




















































with solution 


= C3(A:)(—+ C4(A:)(—+ 


2 (m" - 6) 
H - 4)' 


:E(fc), 

(77) 

where c^{k) and C 4 (A:) are constants. The first two terms 
vanish as t —0 by property P4. This shows that our 
assumption in Eq. ( [6^ is justified. 

From the solutions for the mode functions in Eqs. 
and (771, it is clear that both the adiabatic and the en- 
tropic spectra are proportional to /c^|E(fc)p. Scale in¬ 
variance is obtained if and only if 


which is scale-invariant when X = fi. For a given p. > Vb, 
any ns < (3^^ — 2)/(/r^ — 2) can be achieved by choosing 
A = (ns - l)/fi - n{ns- 3) /2. 

Note that since “just-so” growth implies /(tr) = e~^'^ 
(and property P2 requires/i(tT) = the Lagrangian 

density is given by 


^ = Ir- l^daf - ^e-^"(a</.)2 - (84) 


which is equivalent to the Lagrangian density in Eq. 
with the identifications 


(48) 


|E(fc)p=efc 


-3 


(78) 


for some constant ^ that is independent of k. Eq. (78) is a 
first order differential equation for the coupling function 
/, namely 


4 r o.. 8 

i-r)f (-r)2’ 


its solution is given by 

/(t) = 


sec^ [ln(T/ro)] 
2^(-r)2 ^ 


(79) 


(80) 


where tq is an integration constant. This is clearly not 
monotonic as r —>■ 0 and, therefore, violates property PI. 
Hence, we conclude scale-invariance is impossible for fast¬ 
growing /. 


B. Slow growth: \f,a /f\ < 1 


If I/,<t//| < 1, then 1/7/1 < l/(-r). To leading order, 
both mode functions satisfy 


O 1 

2(^2 _ 2)2 0 ’ 


(81) 


as is clear from Eqs. (52), (53) and ( Mp6 ). In this case, 
both the adiabatic and the entropic spectra are given by 


ns =4- 


y? — Q 
^2-2 


(82) 


which is blue by property P4. Hence, we conclude scale- 
invariance is impossible for slow-growing /. 


C. “Just-so” growth: |/,ct//| 


If l/jo- //I ~ Ij /(o') = 6“'^°' for some A S K such that 
sgn(A) = sgn(/r). The coupling functions a,/3,7,(5 are all 
proportional to l/v7 so the right sides of Eqs. (52) and 


can be neglected. The mode functions effectively 
decouple and the adiabatic spectral index is again given 
by Eq. (82); the entropic spectral index is 


ns = 4- 


Xu — 2 

+ 1 


(83) 


a < —7, 

(j) i —s> X. (85) 

V. DISCUSSION 

In this paper, we have presented explicit examples of 
ekpyrotic models, with and without shift symmetry, that 
have fixed-curve attractor background cosmological so¬ 
lutions generating a scale-invariant spectrum of entropic 
perturbations. The existence of a negative eigenvalue 
associated with the linearized dynamical equations, as 
identified by Tolley and Wesley, can indicate a true insta¬ 
bility requiring fine-tuning of initial conditions for some 
actions. But, for actions of the type discussed here, the 
negative eigenvalue only indicates that the attractor solu¬ 
tion is a curve rather than a point. Fine-tuning of initial 
conditions is thus avoided. Furthermore, as in all ekpy¬ 
rotic models, this class of actions avoids the multiverse 
and the problem that all cosmological outcomes are pos¬ 
sible. Hence, the predictions are “generic,” the same on 
average for any Hubble-sized patch. 

One might be concerned if, in order to avoid fine-tuning 
of initial conditions, the models had to be made more 
complicated. However, the opposite is the case here. The 
actions impose less stringent constraints on the equation- 
of-state, e, during the contracting phase and, hence, less 
fine-tuning of parameters. For such a choice of param¬ 
eters, a scalar spectrum of density perturbations with 
the observed spectral tilt can be generated. As in all 
ekpyrotic models, though, it is not possible to generate 
a detectable spectrum of primordial gravitational waves 
(the ratio of the tensor-perturbation amplitude to the 
scalar-perturbation amplitude, r « 0 ), consistent with 
current limits [misi]. The same class of actions has also 
been shown to generate zero non-Gaussianity during the 
ekpyrotic contraction phase; a small amount of local non- 
Gaussianity may be generated during the bounce, but at 
a level well within current observational bounds on /nl 
[321 [33]. 

It is notable that current cosmological observations are 
constraining ekpyrotic models to be in a class that is the 
simplest, as measured by parameters, degrees of freedom, 
and initial conditions. By contrast, the same observa¬ 
tions are pointing away from the simplest models of in¬ 
flation [51] . 
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VI. APPENDIX 


Here we derive the four properties listed at the begin¬ 
ning of Sec. |IV| using a series of lemmas: 

Lemma 1: a: = 0 cannot generate a scale-invariant adi¬ 
abatic spectrum without fine-tuning. 


li X = 0, z = Zq = const. Since z = a, this means that 
a = zq along the solution. Defining 


0 

b) 

1 

b 

III 

lb 

( 86 ) 

</> = \/ fi^o) h{zo)j , 

(87) 

c 

( 88 ) 

- fid + zo) 

/(=o) ■ 

(89) 

_ h(a + Zo) 

M«) • 

(90) 


the Lagrangian density can be recast as 


£ = ii? - ^(daf - ^F(a)(d^f - VoH(a)e-^'^.(91) 


Eq. (91) has the shift-symmetric form of Eq. ([^ with the 
solution (x, z) = (0, zq) at (x, z) — ( 0 , 0 ) (with z = a) for 
which Tolley and Wesley proved that any solutions of in¬ 
terest are unstable. For fixed-point solutions like this (as 
opposed to fixed-curve solutions) instability implies the 
need for finely-tuned initial conditions. Thus, we con¬ 
clude \a\ oo; we are forced to fixed-curve solutions. 


Lemma 2: / — >■ const cannot generate a scale-invariant 
adiabatic or entropic spectrum without fine-tuning. 


When f = fo = const, the non-canonical coupling be¬ 
comes canonical. Then, the background Eqs. ([^ and ®> 
give {w,x) = (wo,Xo) with xq 7 ^ 0 (cf. Lemma 1). Then 
h = In such canonically-coupled theories, the 

adiabatic perturbation decouples from the entropic per¬ 
turbation and has a blue tilt. More precisely, the equa¬ 
tions of motion to the Lagrangian density 

F = - ^fo(dcf>f - (92) 

linearized around (wo,xo) = (d/\/6, c/\/6fo) yield the 
same equation for both the adiabatic and the entropic 
mode functions 


u 


II 


A 2 ^ 2/0 (c^ + {d^ - 4) fo) 1 \ 

V {c^ + {d^-2)fof (-r)y 


u = 0 (93) 


so that the spectral index is given by 


-k (d^ - 6 ) fo 


2 + (d2 - 2) fo 


(94) 


Avoiding fine-tuning, i.e., negative eigenvalues in the lin¬ 
earized equations of motion, requires / fo d^ > 6 . In 
this regime the spectral indices are blue for both the adi¬ 
abatic and entropic spectra. 

Lemma 3 : If w 0 then f —> {00 or const) at late 
times. 


For scaling solutions, the background equation for w, 
Eq. (§ , can be recast as 

= (95) 


with A = y/6 — l) /x and B = c/w. If A < 0, 

a: < 0 such that z increases with decreasing V; if A > 0, 
a: > 0 and z decreases. Hence, sgn(/') = —sgn(/, 2 ). 
From the right side of Eq. (95), we conclude that, if H < 
0, / —)■ 0 at late times and that, if H > 0, / —)■ (H^ or 00 ) 
at late times. With Lemma 2, we only have to consider 
the case / —>■ 00 . 


Lemma 4: ///—)■ 00 , then w = 0 at late times. 

If / —>■ 00 , the background equation for w in Eq. ([^ 
becomes 


w,N = 3w ( (a:^ -k — I) 


f Vq) 


(96) 


Assume w ^ 0. Then f,z/f = •\/6(a;^ -k — l)/a; and, 
hence, /(cr) = e~^'^ with A = —•\/6(a;^ -\- — l)/x. In 

particular, A and x must have opposite sign. From z,n = 
V6x, we see, though, that x and z must have different 
sign (since N decreases as the system evolves, z will, for 
example, decrease when z,n > 0). Therefore, A and z 
must have the same sign. But then / = e~^^ —7 0 at late 
times, which contradicts our assumption, / —>■ 00 . 

To this point, we have shown that any ekpyrotic, scal¬ 
ing attractor that generates a scale-invariant spectrum of 
either adiabatic or entropic perturbations lies at w = 0 
with / —)• 00 and x 0 . Defining p = y/&x, Eq. ^ 
implies h,^ /h = —/i, from which we conclude, given that 
z = tr, h{a) oc e~^°' (property P 2 ). 

Lemma 5: The scaling solution must correspond to the 
limit given in property PS. 

With /i(cr) = the only scaling solutions of Eqs. ([^ 

and ([^ are {w,x) = {(0, ±1), (0, ^/V6)} if f,a/f 7^ 
const. For the first two solutions, e = 3 . Since ekpy- 
rosis corresponds to e > 3 , we only consider the third 
solution (property P 3 ). 

Note, if f,a /f = const = —A, there is another scaling 
solution at 


{w,x) 


/ , + -6 
1, |A + m| 



(97) 
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At late times, this solution is equivalent to one in Ref. [TS] 
that was shown never to be an attractor. 

Lemma 6 : The coupling /(cr) must grow monotonically. 

Since the solution in property P3 is a hxed-curve, it 
can be parameterized by one variable, 2 . For any fi¬ 
nite \z\, f will be finite so that even if the system lies at 
(w,x) = (0,fi/\/6), the kinetic energy of the ekpyrotic 
field evolves as w,n= — 3 c ( yLt ^/6 — l)v7 (cf. Eq. ([^). 
Thus, in any interval over which / shrinks, \w^n \ grows. 
For this reason, we only consider solutions for which / 


grows monotonically (property P2). Linearizing Eqs. ([^ 
and (|^ about the background in property P3 yields 

<^w,Ar= ^3 - 1^ - (5a;, (98) 

6x,n = 5 - 1 ^ Sx, (99) 

where Sw = w and 6x = x — The eigenvalues 

of this system are {(/r^ — 6)/2, — 6 —/r/,^//)/2}. 

Since fif,z // < 0, this reduced 2x2 system is stable, 
i.e., has positive eigenvalues, if |^| > \/6 (property P4). 
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